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Abstract —Power flow solvable boundary plays an important 
role in contingency analysis, security assessment, and planning 
processes. However, to construct the real solvable boundary in 
multidimensional parameter space is burdensome and time con¬ 
suming. In this paper, we develop a new technique to approximate 
the solvable boundary of distribution systems based on Banach 
fixed point theorem. Not only the new technique is fast and non¬ 
iterative, but also the approximated boundary is more valuable 
to system operators in the sense that it is closer to the feasible 
region. Moreover, a simple solvable criterion is also introduced 
that can serve as a security constraint in various planning and 
operational problems. 

Index Terms —Distribution systems, power flow, solvability, 
security constraint, voltage stability. 

I. Introduction 

The penetration of renewable energy sources is expected 
to increase substantially in the near future. Besides numerous 
advantages of such more environmentally friendly resources, 
additional risks may make the distribution system more vulner¬ 
able to instability fl]. Voltage security is one of main existing 
issues that we need to address carefully and thoroughly in 
order to operate the systems in a reliable and secure fashion. 
Hence, the system operators also need to know the system lim¬ 
its which require faster and more reliable security assessment 
tools. In fact, static voltage stability criteria is widely used by 
system operators 0, 0; moreover, it has been argued that 
static analysis is preferred over dynamic approach 0. Most of 
static analysis techniques rely on the solution of quasi-steady 
state analysis or power flow problem. The concept of quasi¬ 
steady state is based on the assumption of slow dynamics so 
that the power system can be considered in steady state. 

These power flow equations have been posing much diffi¬ 
culty to the engineers due to their inherent nonlinearity and 
high computational cost of finding solutions with conventional 
approaches 0. Iterative methods such as Newton-Raphson 
or Gauss-Seidel methods are the most widely used by power 
flow solvers for the following reasons. First, the solution can 
be found after several iterations in most of the cases, and the 
iterative methods can be easily scaled to bulk systems. Second, 
the solutions obtained can be used for control purposes. How¬ 
ever, these iterative methods may be subject to divergence due 
to either bad initial guesses, proximity to solvable boundaries, 
or the nonexistence of solutions a. 


In order to overcome such handicaps, numerous attempts 
have been made to design more efficient and accurate ways 
to solve or certify existence of solutions to those equations 
0-(9j. Among them, in [10), Bolognani proposed a method 
of certifying the solvability of power flow equations using the 
Banach Fixed Point Theorem. The resulting criterion has a 
simple algebraic representation and dramatically reduces the 
computational of assessing system security. The downside of 
the approach is the conservativeness of the criterion and its 
inability to certify the solvability of the whole region. The 
motivation for this paper is to extend Bolognanis criterion 
and alleviate the conservativeness observed in the original 
definition. 

The key contributions of this paper are as follows. We 
analyze an extension of Bolognani’s criterion and develop a 
new mathematical criterion for solvability certification. Then, 
we derive a simple algebraic representation and propose a 
new technique for fast and non-iterative characterization and 
visualization of the solvability boundaries in a multidimen¬ 
sional parameter space. Finally, two important applications 
for solvable boundary approximation and security constraint 
optimization are presented. 

In section [Til we describe the nonlinear power flow equa¬ 
tions and the notations that be used throughout this paper, 
followed by a brief explanation of the criterion established by 
Bolognani which is the motivation for this paper. The method 
for extending Bolognanis criterion, along with its proof, be 
presented in section m Then, section [IV] displays the 2- 
bus simulation results and related empirical observation that 
discuss the relationship between the resulting ellipsoids and 
rhombuses. This section also discusses the algorithm, or the 
convex hull method for constructing outermost rhombuses that 
can be observed from simulations. Section [V] provides with 
more simulations of 3-bus and 13-bus radial systems, which 
is followed by section m where main results are discussed. 
In section I VIII we introduce two possible application of this 
modification. 

II. Power flow problem and solvable boundary 

APPROXIMATION 
A. Power Flow Equations 

In this work we use the notation introduced originally by 
Bolognani in Ho). Fet C := {l,...,n} be the index set 


denoting all load buses, and let index 0 correspond to the 
slack bus. For h G C, let % £ C and ih £ C be the voltage 
and current, respectively, that define each load bus, and let 
v, i G C n+1 be complex vectors with entries of Vh and ih, 
respectively. Also, let s/, £ C be the imposed power at each 
load bus with s G C” +1 as the vector with Sh as entries. We 
assume that the voltage level at slack bus is fixed at vq = Vq. 

The power flow equations establish the relation between 
the power flows injected and consumed at each bus and the 
corresponding voltage and current levels: 

Vhih = s h Vh G C := 


III. Extension of the solvability certificate 

The solvability certificate defined by (0 is not unique. In 
fact, it was noted in the original work □a that it can be 
generalized by rescaling the power and voltage vectors in the 
original power flow equations. In this work we explore this 
freedom and show that these rescalings can be used to extend 
the criterion, and make it both less conservative and more 
computationally tractable. 

To extend the criterion we replace / with A ip, where A is 
some diagonal matrix in [R" xn and <p some vertical vector in 
[R n . Then we have: 


where the bar denotes the complex conjugate of the given 
vector. The power flow equations need to be supplemented 
with Kirchhoff laws that can be represented in the following 
form: 

vl = vq! + Zip, where Z := 

where 1 is a vector with all components equal to 1, i.e. 1 h = 1. 

In this work, we only consider PQ buses, so the values of 
s’s are assumed to be given. The extension for PV buses is 
rather straightforward and is also described in (10|. 

B. Bolognani’s certificate of solvability 

A sufficient condition for existence of a particular solution 
to the nonlinear power flow equations was derived in Col 
and forms the foundation of the algorithm proposed in this 
work. The criterion is based on reformulation of power flow 
equation as fixed point of nonlinear map, i.e. equation of the 
form / = G(f), where the map is defined via the following 
function: 

G(f) ~ + s c )Z(f+ s c ) (1) 

v o 

and the quantity / is related to bus voltage and current levels 
via the following relation: 


f = Aip = 


\\ipi 

\2<P2 


X n iPn 


It follows from 0 that 


1 \ 

( Pi = ~y2^Pi + \ s i) 'y ' (ZijXj) (fij + Xj Sj) (6) 
0 j =1 

Notice that this equation has the same structure as (0, so 
the same techniques that were used to prove the solvability 
based on 0 can be applied to this equation. Then the revised 
criterion can be rewritten as 


4||ZA|r • IIA-^II < V 0 2 (7) 

The key idea behind the approach explored in this work is 
to characterize the whole union of certificates defined by the 
inequalities 0 for all possible diagonal matrices A. More 
precisely we are interested in characterization of the set 
defined formally as 

V = {s c : 3A such that 4||ZA||* < i||A- 1 s z: || < V 0 2 } (8) 


/ := v 0 i c -sc = V 0 i c ~ s c 

The existence of solution is certified by proving that the 
function G(f) is contractive, and thus the solution exists. The 
final criterion can be compactly expressed as the following 
condition on the power flow levels: 

4||Z||M| S£ || 2 < y 0 2 , (2) 

with the “nuclear” matrix norm defined as 

\\ A \\* = max ( 3 ) 

Similar relation can be derived also for 1-norms. In this case 
one has 

411^11^ • IMIi < V 0 2 (4) 

fl^llSo = maxIAfcl (5) 

Whenever the relation 0 is satisfied for some power flow 
injection vector s, the power flow equations a proven to have 
a solution. Some additional properties of this solution were 
proven in the original work Sla¬ 


in the following sections we present the results of our analysis. 
In the section[IV]we report the results of numerical analysis on 
various models, whereas in section m we provide preliminary 
results from theoretical characterization of the set V. 

IV. Simulation results 

We start our analysis by considering simplest possible radial 
system composed of 2 and 3 buses. To visually characterize the 
union of certificates generated by different matrices A in 0 
we can simply plot the boundaries corresponding to different 
A on top of each other. 

In Figure Q] we present results of analysis of a simple two 
bus system with one slack bus and one PQ bus characterized 
by levels of power consumption defined by P and Q. In this 
model we assume, without loss of generality, that the line is 
purely resistive. More realistic situations can be recovered by 
rotating the vector s in complex plane by an angle defined 
by the impedance of the lines. We overlay the plots of the 
actual solution boundary and the certificates constructed with 
the original criterion © as well as the modified ones 0. 
Several observations can be made from Figure |T| First, it 
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Fig. 1: The dotted curves illustrate the solvable boundary 
from varying rescalings for a 2-bus network. The rescaled 
regions cover only fractions of the solution region certified 
by Bolognani’s criterion, hence represents more conservative 
solution boundaries. 


can be clearly seen that certificates defined by (0 and 0 
do not cover the whole region where the solution to power 
flow equation exists. Second, one can see that rescaling does 
affect the region certified by the criterion. Finally, one can 
see that the original criterion © is the least conservative 
boundary, so rescaling does not increase the size of the region 
certified to have a power flow solution. As we will see, this last 
observation does not apply to more general systems discussed 
below. 

Next, we consider a more complicated case of a 3-bus radial 
system with two PQ buses. In the 3-bus system, both of the 
power lines between buses 1 — 2 and 2 — 3 have the same 
values of resistance and inductance: r = 0.0734p.it. and x = 
0.2581 p.u.. We consider the case in which Pi = P 2 and Q 1 = 
Q 2 - To understand the impact of the rescaling we consider 
diagonal matrices A with elements in the range An £ [0.5, 25], 
The impedance matrix Z in this problem is given by 


Z = 


-I _J L + jJ . 

L 7 203 ^ J 29J 


Figures [2] and Q] illustrate unions of P v’.,v. Q plots with 
regions corresponding to certificates produced by different 
diagonal matrices A. Figure [2] represents the union of 64- 
many PQ plots corresponding to || • H 2 norm as and Figure 
[3] represents the union of plots from || • | [^ norm. The key 
observation that one can draw from these plots is that the 
outer boundary of the union of all certificates has a simple 
polytope, that we refer to as a “Rhombus” in our figures 
due to the shape of its projection. Second, as one can see 
the outer boundary improved the certificate provided by each 



Fig. 2: 64 elliptical solvable regions for the given 3-bus 
network. Note that the contour of the union of the ellipses 
forms a geometry identical to that of the union of rhombuses 
in Figure [3] 
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Fig. 3: 64 rhombus-shaped solvable regions for the 3-bus 
network. 


individual A and provides a less conservative way of checking 
the system solvability. Finally, it is important to note, that the 
outer boundary of both the 11 • 11 J^-norm union and 11 • 11 2 -norm 
union is the same. This observation suggests that this boundary 
can be characterized by a simple algebraic expression. 

On figure [4] we provide yet another empirical confirmation 
of our observation by looking at the projections plots of the 13- 
bus network, assuming P\ = P2 = P 3 = P,\ = P 5 = -y/’o = 
5-P7 = \P& = \Pq = \Pw = j-Pn = \Pi 2 an d Qi = 
Q2 = Q3 = Qa = Qb = ^<36 = \Qi = = 5Q9 = 

10 = jQ 11 = -y Q 1 2 ■ The 13-bus network configuration 
and data is described in CD. Moreover, fixed shunt capacitors 
are installed at buses 3, 4, 6 , 9, 10, 13 with the capacities of 
0.2, 0.1, 0.05, 0.05, 0.05, 0.1, respectively. The unit of the 
capacities are in p.u.. The impedance matrix Z £ |R 12x12 is 
then computed, and the plot for P v.s. Q is plotted using the 
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Fig. 4: Projection of 1 st and 2 nd variables’ plane for 13-bus 
network. On the left is the rescaled regions from the definition 
from ( 0 , and on the right is from 0 . 


modified criterion, using both norms || • ||£ and || • H^. Again, 
we see that the union is always a simple geometric shape that 
looks like a rhombus on the projection plane. 

In the following section, based on our empirical observa¬ 
tions we introduce a conjecture that defines the solvability 
certificate in terms of of the simple piecewise linear function 
defining the convex hull of the critical points. 


V. Constructing the outermost rhombus 


In this section, we propose a method for constructing an 
algebraic representation of the “Maximal rhombus” shape that 
we observed in our numerical experiments. This shape is most 
easily constructed using the representation 0. Consider the 
rescaling defined by X^ 1 = max^ \Zhk\- In this case the 
nuclear norm of the matrix HZAH^ = 1 that can be easily 
checked from the definition 0. The certificate is then reduced 
to 


E 


Sk 

„max 

*fc 


< 1 


(9) 


where the value of .sj. nax is given by the following expression: 


„max 

s k 


V 2 

v 0 


4ma Xh \Zhk\ 


( 10 ) 


Note that the expression 0 defines the convex hull of critical 
points = s™ ax efc, where the unit vectors e*. define the 
rectangular coordinate system in the s-space. One can easily 
show that the polytope defined by 0 is the envelope of union 
of polytopes defined by the rescaled version of 0 . It is also 
easy to check that it encompasses all the critical ellipsoids 
defined by the rescaled version of 0 with A& —»• 00. At this 
stage we don’t know the conditions under which the union 
of all the rescaled ellipsoids coincides with the convex hull 
defined by 0 , however our numerical experiments show that 
it is often the case. 


Fig. 5: The plot displays two different solvable boundaries 
of 3-bus system; one, the actual solution boundary for a 3- 
bus network, and two, the contour of the union of elliptical 
rescaled solvable boundaries. The contour was plotted using 
the method for constructing the ’’Outermost Rhombus”. 


VI. Comparison with real boundary 

As the certificate can be naturally used to guarantee the 
existence of power flow equation solutions it is reasonable to 
assess its conservativeness by comparing the certified region 
to the actual solution boundary. On Figure 0 we perform this 
comparison for the 3-bus example system. Here the outer¬ 
most curve indicates the actual solution boundary found via 
Groebner basis approach explained in 0 . The polytope shown 
represents the first quadrant part of the convex hull as defined 
in ©. As expected, the polytope lies inside the solvability 
region, and provides a remarkably tight approximation of this 
boundary. 



Fig. 6 : The PV curves at load bus 2 with different levels of 
reactive power in the 3-bus system. 

As a supplement to Figure[5] Figure [ 6 ] presents 3 PV curves 
corresponding to different fixed values of reactive power Q'i = 
Q 3 . Pa and Pe indicate the actual value and estimated value 
of maximum active power, respectively. For Q 2 = 0, the two 
values match as shown at point B in Figure 0 Interestingly 
enough, the certificate is tight for purely active power flows 





















or DC power flow, but is more conservative as the amount of 
reactive power consumption is increased. 

VII. Practical applications 

In this section, we propose two important applications of 
the solvability certificate derived above. 

A. Non-iterative approximation of solvable boundaries 

In steady state analysis or power flow problem, the solvable 
boundary consists points corresponding to saddle-node bifur¬ 
cation where the Jacobian matrix becomes singular. Knowing 
the solvable boundary of the system is essential to system 
operators because it indicates the limit of the system and 
the threats to viability of the system operation ca. How¬ 
ever, constructing the solvable boundary in multi-dimensional 
parameter space is difficult due to the complexity of the 
system fl3l . It has been drawing much attention for decades, 
and most of the proposed techniques are computationally 
prohibitive iterative methods lfl~2l . H 1 41 — [fl~6i . which may suffer 
from the divergence problem. On the other hand, the simple 
linear criterion derive in this work a reasonable alternative 
to characterize non-iteratively the solvable boundary. Besides 
computational efficiency provided, a prompt technique for 
visualization of the solvable boundary is extremely helpful 
not only in contingency analysis and planning processes, but 
also in emergency controls where faster decision-making is 
required. Moreover, by construction, the approximated region 
bounded by the outermost rhombus is inscribed in the actual 
solvability boundary, hence, the approximated boundary is 
closer to the feasible region in which the operational con¬ 
straints are satisfied. 

B. Solvability criterion for security constraints 

Another important application can be developed from sim¬ 
ple the solvability criterion described in ©. The simple 
criterion can be incorporated as a security constraint in plan¬ 
ning and operational problems such as OPF. Whenever the 
solvability criterion is satisfied, the system is guaranteed to 
possess a solution to the corresponding power flow problem. 
The simple criterion ([9]) is also useful for corrective control. 
For example, one can use this to find optimal load shedding. 

Moreover, note that the proposed simple criterion can be 
expressed in terms of the slack bus voltage, Vo. Other forms of 
solvable criteria are described in IfTTl . fT8l . In the future work, 
we plan to extend the proposed criterion to take in account 
the feasibility of the solutions, including voltage and current 
criteria. We expect that the new criterion could be applied 
directly to assess quasi-static voltage stability and would help 
to avoid voltage collapse in static scenarios. 

VIII. Conclusions 

In this work, we extended the solvability criterion derived 
in Col by exploring the freedom associated with rescaling 
matrices A. We observed, that the union of the regions 
described by these certificate has often a simple polynomial 
shape and can be represented as a convex hull of relatively 


small number of points represented by the expression (|9|. 
This simple solvability certificate depends only on the system 
configuration and can be naturally used as voltage security 
constraint. 

In future works, we plan to develop a new method to 
identify feasible load change and design secondary voltage 
control based on the definition of the polytopic solvability 
certificate. Moreover, we plan to explore the possibility of 
using this criterion as static voltage stability indicator. 
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